Hadronic resonance gas and charged particle's pt spectra and elliptic flow in y^=200 

GeV Au+Au collisions 
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Charged particles pr spectra and elliptic flow in 0-60% Au+Au collisions at RHIC are analyzed in 
a hydrodynamic model with hot, hadronic resonance gas in the initial state. Physically conceivable 
hadronic resonance gas, thermalized in the time scale Ti=l fm, at a (central) temperature Ti— 220 
MeV, with viscosity to entropy ratio ri/s— 0.24, reasonably weU explains the pr spectra in all the 
collision centralities. Centrality dependence of elliptic flow however demands continual increase of 
viscosity to entropy ratio with centrality. 



Recent experiments in Au+Au collisions at Relativis- 
tic Heavy Ion Collider (RHIC) [1],;2J,_3J,[4] produced con- 
vincing evidences that in central and mid central Au+Au 
collisions, a hot, dense, strongly interacting matter is cre- 
ated. Whether the matter can be identified as the lattice 
QCD [5] predicted Quark Gluon Plasma(QGP) or not is 
still a question of debate. Confirmatory identification of 
the matter with QGP requires complete elimination of 
the possibility that a hadronic state is produced in the 
initial collisions. It is yet to be done. Information about 
the initial state is always indirect. QGP is a transient 
state and even if produced, exists for a short time scale, 
it expands, cools and eventually transforms into hadrons. 
Hadrons are the experimental observables and any infor- 
mation about the initial state has to be obtained from 
the observed hadrons. Dynamical models are essential 
to extract information about the initial state from the 
experimental observables. 

Relativistic hydrodynamics is one of the few dynami- 
cal models which has been successfully applied in RHIC 
collisions to obtain information about the initial medium 
produced in Au+Au collisions. It is assumed that in a 
Au+Au collision a fireball is produced. Constituents of 
the fireball collide frequently to establish local thermal 
equilibrium sufficiently fast and after a certain time r^, 
hydrodynamics become applicable. If the macroscopic 
properties of the fiuid e.g. energy density, pressure, ve- 
locity etc. are known at the equilibration time Ti, the 
relativistic hydrodynamic equations can be solved to give 
the space-time evolution of the fireball till a given freeze- 
out condition such that interactions between the con- 
stituents become too weak to continue the evolution. Us- 
ing suitable algorithm (e.g. Cooper-Frye) information at 
the freeze-out can be converted into particle spectra and 
can be directly compared with experimental data. Thus, 
hydrodynamics, in an indirect way, can characterize the 
initial state of the medium produced in heavy ion colli- 
sions. Hydrodynamic equations are closed only with an 
equation of state and one can investigate the possibility 



of phase transition in the medium. 

Relativistic hydrodynamics with QGP in the initial 
state, thermalized in the time scale tj wO.6 fm, at an 
initial central energy density en « 30 GeV/ Jm^, could 
explain a host of experimental data produced in y^=200 
GeV Au+Au collisions 6]. In an alternative to the ini- 
tial QGP state, namely the hot hadronic resonance gas 
(HRG), description to the experimental data gets much 
poorer 7]. Also, to reproduce the experimental multiplic- 
ity, hot hadronic resonance gas is required to be initial- 
ized at very high temperature, e.g. Ti « 270 MeV. Den- 
sity of HRG grows rapidly with temperature. At T « 270 
MeV, density of hadrons is large, phad ~ 4 fm~^. At such 
a large density hadrons will overlap extensively and it is 
difficult to believe that they could retain their individual 
identity. It does appear that a consistent description of 
RHIC data could not be obtained in a hadronic resonance 
gas model. However, all the analysis of RHIC data with 
hadronic resonance gas in the initial state was performed 
in the ideal hydrodynamic limit. In recent years, there is 
much proCTess in pra ctical application of viscous hydro- 
dynamics B, iL UMIi, M: MM, M, 0: [13, 
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Several codes were developed to 
numerically solve causal dissipative hydrodynamics and 
have been applied successfully to analyze experimental 
data at RHIC collisions with QGP in the initial state. 
However, we have not come across any analysis of RHIC 
data with 'viscous' hadronic resonance gas in the ini- 
tial state. Viscous effects can be large in hadronic reso- 
nance gas. Model calculations t2J| , [2J| , [2g] , [2g] indicate 
that viscosity to entropy ratio of a hadronic resonance 
gas can be considerably larger than the ADS/CFT limit, 
v/s > l/47r ^3l- Since entropy is generated during evo- 
lution, unlike an 'ideal' HRG, a 'viscous' HRG can be 
initialized at a lower temperature such that hadrons re- 
tain their identity yet reproduce the experimental mul- 
tiplicity. Indeed, in one of the earliest applications of 
one dimensional viscous hydrodynamics, it was shown 
that WA80 single photon spectra, though over predicted 
in an ideal hadron gas evolution, is reasonably well ex- 
plained when viscous effects are accounted for |23|- To 
exclude the possibility of hadron gas formation in initial 
Au+Au collisions, instead of QGP, it is important that 
RHIC data are analyzed in a viscous hadronic resonance 



gas model. 

We assume that in Au+Au collisions a baryon free, 
hot hadronic resonance gas, comprising all the hadron 
resonances with mass rrires <2.5 GeV is produced. The 
hot, hadronic resonance gas is assumed to thermalize in 
the time scale Ti~l fm. We also neglect all the dissipative 
effects except for the shear viscosity. The space time 
evolution of the fluid is obtained by solving. 



gas. For hadron size « 0.5 fm, limiting hadron density 
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Eq[T] is the conservation equation for the energy- 
momentum tensor, T^"" = (e + p)u''u'' - pg^"' + vr^'', 
e, p and u being the energy density, pressure and fluid 
velocity respectively, tt^'^ is the shear stress tensor (we 
have neglected bulk viscosity and heat conduction) . Eql2] 
is the relaxation equation for the shear stress tensor tt^"^ . 
In Eql21 D = u^df^i is the convective time derivative, 
y<t^u''> = \{Wu'' + Vu^) - \{^.u){g^''' - ut^u") is 
a symmetric traceless tensor, rj is the shear viscosity 
and Ttt is the relaxation time. It may be mentioned 
that in a conformally symmetric fluid relaxation equation 
can contain additional terms }17l] . Assuming longitudinal 
boost-invariance, the equations are solved with the code 
'AZHYDRO-KOLKATA'fUl in (t = Vt^ - z^,x,y,r] = 
■^Injz^) coordinates. 

Solution of EqU] and [2] require initial energy density, 
velocity distribution in the transverse plane at the ini- 
tial time. We assume that at the initial time Tj=1 fm, 
initial fluid velocity is zero, Vx{x,y) = Vy{x,y)=0. The 
initial energy density in an impact parameter b collision 
is assumed to be distributed as fal. 



£0^,x,y) ^ eo[{l- x)Npart{^,x,y) + xNcou{h,x,y)] (3) 

where NpartO^jXjy) and Ncou{h,x,y) are the transverse 
profile of participant numbers and binary collision num- 
bers respectively. Nparti^), x, y) and Ncoiii)^, x, y) can be 
calculated in a Glauber model, x in Eql3]is the fraction of 
hard scattering. Most of the hydrodynamic simulations 
are performed with hard scattering fraction a;=0.25 or 
Recently, in [32| 



0.13 



31 



it was shown that in colli- 
sions beyond 0-10% centrality, simultaneous description 
of charged particles px spectra and elliptic flow are best 
obtained with hard scattering fraction x — Q. Only in 0- 
10% centrality collisions, hard scattering fraction a; = 1 is 
preferred. In the following, in all the collision centralities, 
we assume hard scattering fraction x—Q^ understanding 
that we may under predict elliptic flow in very central 
collisions. Eq in Eql2]is the central energy density of the 
fluid in an impact parameter b=0 collision. Generally, eg 
is obtained by fitting experimental data, e.g. multiplicity, 
Pt spectra etc. We however fix Eq to the highest possible 
value for a 'physically' conceivable hadronic resonance 



(such that hadrons are not overlapped) is Pi^^t = l/V 



2 fm '*. pf°:^.^—2 fm "^ corresponds to limiting tem- 
perature or energy density, Turmt^ 220 MeV, eiimit=5.1 



GeV/fr 



We fix the central energy density at this 



limiting value, £o = £/imit=5.1 GeV/ fm^. Recent lat- 
tice simulations 'S', with almost physical quark masses 
(ttItt ~220 MeV), predicts a confinement-deconfinement 
transition temperature Tc = 196 ± 3 MeV. Hadrons can 
exists at T w200 MeV. Existence of HRG at -10% higher 
temperature r=220MeV is a definite possibility. 

Dissipative hydrodynamics also require initialization of 
the shear stress tensor tt'^" as well as the relaxation time 
Tjr. We assume that shear stress tensors are initialized 
at the boost-invariant values, tt^^ = tt^^ = 2r]{x,y)/3Ti, 
TT^^ = [1^]. For the relaxation time Ttt, we use the 
Boltzmann estimate, t^ = Gij/Ap w ^^. Finally, hy- 
drodynamic models require a freeze-out condition. We 
assume that the fluid freeze-out at a fixed temperature 
Tf=110 MeV. 

Only unspecified parameter in the model is the viscos- 
ity coefficient r]. We assume that throughout the evolu- 
tion, the ratio of viscosity to entropy density rj/s remain 
a constant and simulate Au-|-Au collisions for four val- 
ues of ri/s, (i) T]/s=0 (ideal fluid), (ii) ry/s = l/47r=0.08 
(ADS/CFT lower limit of viscosity) and (iii) r//s=0.16 
and (iv) ?7/s=0.24. Assumption of constant rj/s flxes the 
variation of viscosity with temperature, yy ex s oc T^. 
In FigH] and [U simulated charged particles pr spec- 
tra and elliptic flow are compared with experimental 
data. In FigUl in six panels, PHENIX data ^29] for 
charged particles pr spectra in 0-10%, 10-20%, 20-30%, 
30-40%, 40-50% and 50-60% Au-hAu colhsions are shown. 
The dashed, dashed-dot, short dashed and solid lines in 
the figure are simulated spectra from evolution of hot 
hadronic resonance gas, with viscosity to entropy ratio 
ij/s^O, 0.08, 0.16 and 0.24 respectively. If viscous ef- 
fects are neglected, hot hadronic resonance gas, initial- 
ized with central temperature Ti— 220 MeV, do not ex- 
plains the data, the data are largely under predicted. For 
example, at px ~1.75 GeV, in all the collision centrali- 
ties, ideal hot hadronic gas under predict PHENIX data 
by a factor of —6. Data at higher px are even more 
under predicted. When viscous effects are included, par- 
ticle yield increases, more at high px than at low px and 
discrepancy with experiment and simulated spectra di- 
minishes. At Px «1.75 GeV, data are under predicted, 
by a factor of — 4 in evolution of minimally viscous HRG 
(?7/s=0.08) and by a factor of — 2 in evolution of HRG 
with 7y/s=0.16. Experimental data are reasonably repro- 
duced with viscosity to entropy ratio 77/3=0. 24. Indeed, 
in all the centrality ranges of collisions, simulated spectra 
from evolution of hadronic resonance gas with viscosity to 
entropy ratio 77/5=0. 24, agree with experimental within 
—10%. It is very interesting to note that the 77/5=0.24, 
obtained from the analysis is in close agreement with the- 
oretical estimates of viscosity to entro py ratio of a hot 
hadronic resonance gas, 77/5=0. 24-0. 30 [23l.l26|. 
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FIG. 1: In six panels, PHENIX measurements 29] for charged 
particles pr spectra in 0-10%, 10-20%, 20-30%, 30-40%, 40- 
50% and 50-60% Au-f Au collisions are shown. The dashed, 
dash-dotted, short dashed and solid lines are simulated spec- 
tra from evolution hot hadronic resonance gas with rj/s—O, 
r;/s=0.08 and r;/s=0.16 and 77/s=0.24 respectively. 
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FIG. 2: In six panels, PHENIX measurements [30] for 
charged particles elliptic flow in 0-10%, 10-20%, 20-30%, 30- 
40%, 40-50% and 50-60% Au-|-Au collisions are shown. The 
dashed, dash-dotted, short dashed and solid lines are sim- 
ulated flow from evolution hot hadronic resonance gas with 
rj/s—O, rj/s— 0.08 and r;/s=0.16 and ri/s—0.24 respectively. 



The analysis clearly indicate that at least for the 
charged particles pr spectra , it is not necessary that 
QGP fluid is produced in Vs=200 GeV Au-^Au colli- 
sions. 'Physically conceivable' hot hadronic resonance 
gas, with viscosity to entropy ratio 77/5=0. 24, could re- 
produce the spectra. Elliptic flow analysis on the other 
hand gives a different result. In figlH PHENIX measure- 
ments Kid\ for charged particles elliptic flow in 0-10%, 
10-20%', 20-30%, 30-40%, 40-50% and 50-60% Au-^Au 
collisions are shown. PHENIX collaboration measured 
charged particles V2 upto pT ~ 8 GeV. In Figl2l mea- 
surements upto pt—S GeV are shown only. Hydrody- 
namic models are not well suited for high pT particles. 
In order to study non-flow effects that are not correlated 
with the reaction plane, as well as fluctuations of V2, 
PHENIX collaboration obtained V2 from two indepen- 
dent analysis, (i) event plane method from two indepen- 
dent subdetectors, V2{BBC} and V2{ZDC-SMD} and 
(ii) two particle cumulant V2{2}. V2{2} from two parti- 
cle cumulant and V2{BBC} or V2{ZDC — BBC} from 
event plane methods agree within the systematic error. 
It may also be mentioned here that V2{2} in PHENIX is 



lower than that U2{2}in STAR measurements, but they 
agree within the systematic error. All the three measure- 
ments of V2 are shown in Figl^l As before, the dashed, 
dashed-dotted, short-dashed and solid lines in Figl5] are 
simulated flow with ry/s=0, 0.08, 0.16 and 0.24 respec- 
tively. Unlike the px spectra, centrality dependence of 
elliptic flow are not explained with a single value for vis- 
cosity to entropy ratio. Data demand more viscous fluid 
in more peripheral collisions. For example, in 0-10% cen- 
trality collisions, experimental flow are under predicted 
in ideal HRG. With viscous HRG, flow is even more under 
predicted. We have assumed participant scaling for the 
initial energy density. As mentioned earlier, very central 
collisions prefer binary collision number scaling of initial 
energy density 32] . Spatial eccentricity is comparatively 
large in binary collision scaling and elliptic flow in 0-10% 
collision will be better explained if initial energy density 
scales with binary collision numbers. In 10-20% central- 
ity collisions, experimental flow is under predicted with 
77/5=0.16 and 0.24. Data demand HRG with viscosity 
in the range ?7/s=0-0.08. Elliptic flow in 20-30% cen- 
trality collisions prefers more viscous HRG, rj/s ss0.16. 
Hadronic resonance gas with 77/s w 0.24 approximately 



explains flow in 30-40% centrality collisions. Elliptic flow 
in 40-50% and 50-60% Au-hAu collisions demand HRG 
with r]/s > 0.24. Continual increase of viscosity with cen- 
trality can be understood. In a hadronic resonance gas 
viscosity to entropy ratio increases with decreasing tem- 
perature. Initial temperature of the fluid also decreases 
as the collisions become more and more peripheral. For 
example, central temperature of the fluid is Ti w 219 
MeV in a 0-10% collision. In a 30-40% collision, cen- 
tral temperature is Ti «210 MeV. Hadronic fluid will be 
more viscous in a peripheral collision than in a central 
collision. Elliptic flow, being a sensitive observable can 
detect the small variation in temperature. 

It is however difficult to explain why the stated change 
in temperature is not reflected in the px spectra. Though 
centrality dependence of elliptic flow require continual 
increase of viscosity, charged particles px spectra def- 
initely do not demand such increase . Rather the px 
spectra are well explained with a fixed value for viscos- 
ity, r]/s Ri0.24(see FigH]). It appear that a viscous hy- 
drodynamic model, with hadronic resonance gas in the 
initial state, is incapable of simultaneous explanation of 
centrality dependence of charged particles px spectra and 
elliptic flow. We may note here that hydrodynamics, with 
QGP in the initial state, also require continual increase 
of viscosity to explain the centrality dependence of ellip- 
tic flow [2DJ, particle spectra are however are explained 
with a fixed viscosity to entropy 22] . 

Before we summarise, we note that we have used some 
specific initial conditions for the hydrodynamical model, 
e.g. initial (central) energy density Si— 5.1 GeV/fm^, 
initial or thermalization time Ti=l fm, freeze-out tem- 



perature Tf=110 MeV. The initial energy density was 
fixed from physical requirement that constituents of the 
hadronic resonance gas do not overlap extensively and 
loss their identity. Hydrodynamic model analysis of 
RHIC data with QGP in the initial state indicate ther- 
malization time for QGP Ti— 0.6 fm. Compared to QGP, 
a hadronic resonance gas is expected to take longer time 
to achieve thermalization and the thermalization time 
for the hadronic resonance gas was fixed to the canonical 
value Ti=l fm. Similarly, freeze-out temperature Tp=110 
MeV was indicated in the hydrodynamic model analysis 
of RHIC data [6|. All physically supported parameters 
set were not examined. 

In summary, we have studied the possibility of explain- 
ing charged particle's px spectra and elliptic fiow with- 
out invoking Quark-Gluon-Plasma. We assume that a 
hot hadronic resonance gas is formed in initial Au+Au 
collisions. The hot hadronic gas thermalizes in the time 
scale Ti sa 1 fm to central temperature Ti w 220 MeV. 
While an ideal hadronic resonance gas do not explain 
the charged particle px spectra (data are largely un- 
der predicted), a viscous hadronic gas, with viscosity 
to entropy ratio rj/s ssO.24, explains the px spectra in 
all the centrality ranges of collisions. Elliptic flow data 
however are more sensitive and demand more viscous 
HRG in peripheral collisions than in central collisions. 
We conclude that with the initialization and freeze-out 
parametrization presented here, a viscous hydrodynamic 
model containing only hadronic resonance gas, is inca- 
pable of simultaneous explanation of centrality depen- 
dence of charged particles px spectra and elliptic flow in 
Vi=200 GeV Au-|-Au coUisions at RHIC. 
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